Local-density approximation for confined bosons in an optical lattice 
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We investigate local and global properties of the one-dimensional Bose-Hubbard model with an 
external confining potential, describing an atomic condensate in an optical lattice. Using quan- 
tum Monte Carlo techniques we demonstrate that a local-density approximation, which relates the 
unconfined and the confined model, yields quantitatively correct results in most of the interesting 
parameter range. We also examine claims of universal behavior in the confined system, and demon- 
strate the origin of a previously calculated fine structure in the experimentally accessible momentum 
distribution. 
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I. INTRODUCTION 

When a system of cold atoms confined in an external 
trap is exposed to a standing-wave laser field the electric 
field couples to the dipole moments of the atoms and a 
so-called optical lattice it created. A system of trapped 
bosonic atoms in such a lattice is expected to be well de- 
scribed by a Bose-Hubbard models The phase diagram 
of the homogeneous Bose-Hubbard model was obtained 
by Fisher et al.£ and confirmed by quantum Monte Carlo 
investigations,- The model exhibits two possible phases 
for large on-site interaction between the atoms. At com- 
mensurate filling there is a Mott insulating phase and 
at non-commensurate fillings there is a superfluid phase. 
The optical lattices constitute a new realization of the 
Hubbard model, and offer an unprecedented control of 
most model parameters such as lattice size, dimension, 
chemical potential, and the intersite coupling. Using this 
high degree of control Greiner et al. have recently studied 
the release of atoms from an optical latticed In an ab- 
sorption measurement they studied the momentum dis- 
tribution and detected a gradual decrease of the center 
peak in the absorption spectra with an increase in the 
on-site interaction between atoms, indicating a growing 
Mott insulating phase in the center of the trap. The 
need to interpret and understand this, and related, recent 
experiments performed in optical lattices have renewed 
interest in understanding various extensions of the well- 
studied bosonic Hubbard model. 

In the case of optical lattices the picture is compli- 
cated by the trapping potential. The trapping poten- 
tial is an external potential applied to keep the atoms in 
the lattice. Due to the spatial variation of the confin- 
ing potential, different phases can be realized in differ- 
ent parts of the trap, leading to a spatial phase separa- 
tion. This has been demonstrated using quantum Monte 
Carlo techniques for bosons^ and fermions.- 6 ^ These stud- 
ies focused on the differences between the phase dia- 
gram of the confined and unconfined model and calcu- 
lated state-diagrams for the confined model. A density 
matrix renormalization group study^ shows that the de- 
cay of correlation functions in the confined case can be 



easily rescaled to obtain the form expected for the homo- 
geneous model. These results were also shown to agree 
very well with a hydrodynamical treatment of the ID 
Bose gas combined with a local density approach. Results 
for large systems have also been obtained by combining 
the Gutzwillcr mean-field ansatz with a numerical renor- 
malization group procedure^ Recent large-scale Monte 
Carlo studiesifiiii extend previous work to higher dimen- 
sions and address questions of critical behavior and local 
order parameters. 

In this work we focus on the similarities of the phase 
diagram for the unconfined and the confined state dia- 
gram. We map the site-dependent confining potential 
for the trapped system to the chemical potential for the 
homogeneous system and in this manner determine the 
properties of the confined system in a local-density ap- 
proximation. It is reasonable to believe that a good 
agreement between the confined and the unconfined sys- 
tem is obtained for local observables in weakly interacting 
systems, where the correlation length is short. Here we 
demonstrate that the approximation can be used with 
good accuracy in most of the range of interactions where 
the Mott insulating phases appear. Furthermore, the ap- 
proximation can be used also for the nonlocal momen- 
tum distribution, which is calculated from the particle- 
particle correlation function. This local-density approx- 
imation has been used to compare a number of observ- 
ables previously^SiLSiiSiii but here we present a more 
exhaustive comparison of several observables. One prac- 
tical aspect of the local-density approximation is that it 
allows a very quick estimate of the distribution of differ- 
ent phases in a specific trap setting from knowledge of 
only the phase diagram of the homogeneous model. 

The outline of the paper is as follows. In Sec.[H]we dis- 
cuss the phase diagram for the Bose-Hubbard model and 
explain the mapping between the confined and the un- 
confined model. In Sec. IIIII we briefly discuss the Monte 
Carlo method and present results testing the validity of 
the local-density approximation. We conclude with a 
summary and conclusions. 
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II. BOSONS IN AN OPTICAL LATTICE 

Starting from a general Hamilton operator for bosonic 
atoms in an external trapping potential and optical lat- 
tice, Jaksch et alX expand the bosonic field operators 
in the Wannier basis and arrive at the Bose-Hubbard 
model. In order to review the properties of this funda- 
mental model we start by disregarding the effects of the 
confining potential and consider the homogeneous one- 
dimensional Bose-Hubbard model, 

. , , v 

H BH = -t(c\c i+1 + H.c.) + —n^rii - I) - ^n 2 , (1) 

i 

where t is the hopping matrix element between adjacent 
sites, V is the on-site repulsion, fj, is the chemical poten- 
tial, and rii is the number of bosons at site i. At zero tem- 
perature there are two possible phases for the model. For 
small on-site repulsion the kinetic energy is minimized 
by a phase coherent superfluid phase, characterized by a 
nonzero superfluid density, algebraically decaying corre- 
lation functions and a linear, gapless excitation spectrum. 
As the on-site repulsion is increased it becomes harder for 
the bosons to move, due to the increasing potential en- 
ergy cost. At a critical value of the on-site repulsion the 
model enters the Mott insulating phase where each lat- 
tice site is filled with the same number of bosons. The 
suppression of the fluctuations in particle number makes 
the Mott insulating phase incompressible with a vanish- 
ing compressibility, k = f3((n 2 ) — (n) 2 ) = 0. The insu- 
lating phase is gapped and correlation functions decay 
exponentially. The zero-temperature phase diagram for 
the Bose-Hubbard model, as a function of density and 
the hopping parameter t/V, is schematically shown in 
the right part of Fig. ^ For small values of t/V there is 
a series of Mott lobes as the density is increased. In each 
Mott lobe the system is in the insulating phase and the 
density is constant. At large values of t/V the model is 
in the superfluid phased 

As shown by Jaksch et alX the Bose-Hubbard model 
is an effective model for a Bosc-Einstein condensate of 
atoms when subjected to an optical lattice potential. The 
potential energy is given by V = 4ira s h J d 3 x\w(x)\ 4: /m, 
where a s is the s-wave scattering length, w(x) is a local- 
ized Wannier function and m is the atomic mass. As the 
lattice potential is turned on the condensate remains in 
a superfluid state for small values of the potential. How- 
ever, as the potential is increased there is a transition 
to an insulating phase where the phase coherence of the 
condensate is destroyed, as demonstrated by Greiner et 
al.X Other properties of the Mott insulating phase, such 
as the robustness to external perturbations due to the ex- 
citation gap have also been verified experimentally. 4 The 
relevance of the Bose-Hubbard model to Bose-Einstein 
condensates have therefore been demonstrated both ex- 
perimentally and theoretically, and in this work we exam- 
ine the model further using quantum Monte Carlo sim- 
ulations. Some of the quantities we examine, such as 
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FIG. 1: Demonstration of the local-density approximation 
where a system with a confining potential is mapped to a 
slice in the phase diagram for the Bose-Hubbard model. The 
figure to the left represents the site dependent potential /Xi /V 
as a function of lattice site for a confined system with hop- 
ping parameter t/V — t' . The figure to the right shows a 
schematic phase diagram for the unconfined Bose-Hubbard 
model, where the insulating Mott lobes (MI) and superfluid 
phase (SF) are indicated. The dashed line in the phase dia- 
gram represents the different values of the effective chemical 
potential obtained in the confined system. The three x's mark 
how three specific sites in the confined system are mapped to 
the phase diagram on the right. 



the density profile of bosons in the trap, or the momen- 
tum distribution of the bosons, can be directly measured 
experimentally. Other quantities, such as local-density 
fluctuations may be harder to observe experimentally, 
but are of interest since they tell us how the superfluid- 
insulator transition occurs in a confined system. 

In experimental realizations of the optical lattice a con- 
fining potential is necessary to prevent the bosons from 
escaping. The confinement typically spans a couple of 
hundred lattice sites and its strength is adjusted to be 
large enough so that there are no atoms at its edges. The 
confining potential is also used to control the dimension- 
ality of the system. By using a very narrow confinement 
in some directions the dimensionality of the system is 
changed. In this work we consider a setup with a narrow 
confinement in two dimensions and a broad confining po- 
tential in one dimension, leaving a one-dimensional sys- 
tem. The broad confining potential adds a site dependent 
potential term to the Hamiltonian 



H = H 



BH 



N 
~2 



(2) 



V 

-t(c}a +1 + H.c.) + —ni(ni - I) 



where a is a parameter that controls the power of the 
confining potential. In the second line of the equation 
the confining potential and the chemical potential fi are 
combined into a site dependent chemical potential, 
Different sites in the lattice therefore have different effec- 
tive chemical potentials. In the limit of vanishing t the 
local properties of the system are determined by the value 
of the local chemical potential. For finite values of t/V 
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we assume that the local properties of the confined sys- 
tem may be obtained from a homogeneous system with 
the same value of t/V as in the confined system, and a 
chemical potential equal to the effective chemical poten- 
tial at the point considered in the confined lattice. This 
mapping is illustrated in Fig. ^ where it is shown how 
the confined system represents a constant t/V slice in 
the fi — t diagram of the homogeneous system. 

We examine the accuracy of this local-density approxi- 
mation for a range of hopping parameters where the Mott 
insulating phases are present. In general we expect poor 
results for very deep and narrow traps where the poten- 
tial varies rapidly from site to site and close to phase 
transitions where the correlation length diverges in the 
unconfined system. At a true phase transition the cor- 
relation length diverges, but due to the limited size of 
the confined system the correlation length cannot diverge 
and there cannot be a true phase transition. We merely 
expect a system with spatially separated regions of differ- 
ent phases, a scenario which is supported by our results. 
A careful analysis of the lack of critical behavior can also 
be found in Refs. 

The trap used in the experiment described by Greiner 
was of size 65 x 65 x 65 and contained a maximum of 2.5 
bosons per site. In this work we consider one-dimensional 
traps with a linear size of about 300 sites and a density 
up to two bosons per site. Therefore, we believe that the 
results we obtain are pertinent to recent experimental 
setups. 



III. RESULTS 

In this work we use the stochastic series expansion 
quantum Monte Carlo methodi^S*^ We perform the cal- 
culations at a fairly low temperature of (3/V = 256, and 
the convergence is controlled with a calculation at a tem- 
perature of half this value. 

We study the densities in a trap with a linear, square, 
and a quartic potential, i.e. a = 1, 2, and 4 in Eq. 
The calculation is performed for three different values of 
the hopping parameter, t — 0.05, 0.1, and 0.2. We mainly 
show results for the two cases t — 0.1 and 0.2, since t = 
0.05 and t = 0.1 yield quantitatively very similar results. 
The on-site repulsion V is set to one. The strength of 
the confining potential U c , and the chemical potential 
\i, are adjusted so that the effective chemical potential 
Hi/V = 1 in the center of the well, and is large enough to 
ensure that there are no bosons at the edges of the trap. 

In Fig. |3 we display the approximate phase diagram of 
the homogeneous Bose-Hubbard model, based on Fig. 5 
in Ref. and Fig. 9 in Ref. 0] . In the phase diagram 
we indicate the three slices at fixed value of the hopping 
t/V that our confined models represent. We note that 
the slices at t/V = 0.05 and t/V = 0.1 cut through a 
significant part of the first Mott lobe, while t/V — 0.2 is 
chosen to cut through the tip of the same lobe. 

The local-density as a function of lattice site is shown 
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FIG. 2: A detailed phase diagram showing the first Mott 
lobe (MI) and the surrounding superfluid phase (SF). The 
chemical potential fi/V is plotted as a function of the hopping 
parameter t/V . The three slices of the fj, — t diagram that our 
confined systems represent are indicated with the dotted lines. 
The dotted lines start at a p low enough to ensure that there 
are no atoms at the edges of the trap, and they all end at 
H/V = 1. 



in Fig. |2| We note that for £ = 0.1 there is a range of 
lattice sites for which the local-density takes the value 
of one boson per site for all three confining potentials. 
Since the density fluctuations are strongly suppressed in 
these regions they appear incompressible and are there- 
fore called Mott plateaus, indicating that the system is in 
a Mott insulating state. For t — 0.2 there is not yet a well 
defined plateau, but all three density curves correspond- 
ing to the three different confining potentials display a 
slight kink at unit density. 

If the densities at different lattice sites are presented 
as a function of the effective on-site chemical potential, 
this density can be directly compared with the density in 
a homogeneous system with the same value of the chem- 
ical potential. The results are shown in Fig. 01 and we 
see a good quantitative agreement between the densities 
for the confined and the unconfined systems for most 
of the parameter regime. Due to the diverging correla- 
tion function there is some discrepancy close to the place 
where there is a phase separation between the Mott insu- 
lating regions and the superfluid phases, which is shown 
in the insets. The sharp edge for the homogeneous sys- 
tem is simply smoothed out, indicating the absence of a 
true phase transition as the plateau is approached in the 
confined system. In the lower part of the figure, which 
presents results for t/V — 0.2, one can see that there is a 
small plateau in the homogeneous case, which is smeared 
out in the confined systems. This indicates that there 
is no true Mott insulating region in the confined system 
even though the unconfined system has a Mott insulating 
phase. We also can detect a small difference close to the 
edge of the trap, where the local-density decreases quite 
rapidly from site to site. 

The Mott insulating phase is incompressible, and the 
global compressibility vanishes in the Mott insulating 
phase for the unconfined system. For the confined system 
the global compressibility never vanishes* since there are 
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FIG. 3: The local-density m as function of lattice sites for 
three different confining potentials. The upper figure is cal- 
culated for t/V = 0.1 and the lower for t/V = 0.2. The lattice 
constant is denoted by a. 

always superfluid regions present. Here we examine the 
behavior of a local compressibility, defined as the fluctu- 
ations in the local-density 

A? = (n?> - (n,) 2 . (3) 

The results are presented in Fig. [5] where we once again 
compare the result for the confined system with values 
obtained for the unconfined case. Again we notice a good 
general agreement, but discrepancies close to the edge of 
the plateaus and particularly at the edge of the trap. The 
quartic potential rises faster than the other potentials at 
the boundary of the trap, and it is also in the data for 
the quartic potential that we see the largest deviations 
from the homogeneous case. For t = 0.2 the inset shows 
the deviations in the region where the homogeneous sys- 
tem forms a Mott plateau. This is another indication 
that the confined system does not manage to form a well 
developed insulator here, as mentioned above in the argu- 
mentation about the density profile. This has previously 
been noted for the case of fermionic systems^ 

Next we focus on the properties of the system as a 
plateau is approached. As a Mott plateau is entered by 
scanning the chemical potential for the homogeneous sys- 
tem there is a second-order phase transition from the su- 
perfluid to the Mott insulating phase. In the confined 




FIG. 4: A comparison of the densities for the confined and 
unconfined case. The local-density m is plotted as a function 
of the effective chemical potential fa/V. Note that the re- 
sults for the quartic, linear, and square potentials essentially 
coincide with each other, and with the results for the uncon- 
fined system. The insets show the region around points where 
there is a phase transition for the homogeneous system. In 
the upper figure t/V = 0.1 and in the lower figure t/V = 0.2. 



system, at least in one dimension, there can be no true 
phase transition since there cannot be a divergent corre- 
lation length at the boundary between the two phases, 
and we simply have a system displaying spatial phase 
separation. A similar conclusion was also reached for 
a two-dimensional systemiA2*iI Nevertheless, there have 
been claims of universal behavior as a Mott plateau is 
approached^ In Ref. @ it was argued that the vari- 
ance in the local-density shows universal behavior as the 
Mott plateau is approached. This behavior was not re- 
produced in a numerical renormalization group study 9 , 
and in Fig. we display the variance in the local-density 
as the Mott lobe is approached for three different val- 
ues of the hopping parameter. By analyzing the data 
in Fig. H3 more closely we find that the local-density ap- 
proximation works very well, in contrast to the fermionic 
casepi The accuracy of the local-density approximation 
probably explains the "universal" behavior seen previ- 
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FIG. 5: Variance in the local-density Aj as a function of the 
effective potential \n/V for different confining potentials and 
also for the unconfined system. In the upper figure t/V = 0.1 
and in the lower figure t/V = 0.2. 

ously. However, as t is varied the curves readily separate 
and there is thus no real universality. 

Another case of universality as the Mott plateau is 
approached in a confined system has been reported^ for 
a different definition of a local compressibility, 

k\ = (nini +j ) - (ni}{n i+j ), (4) 
U\<i 

which reflects the response at site i to a change in the 
chemical potential in a region of size I. This compress- 
ibility was introduced as a local order parameter, which 
vanishes in the insulating phase. It was also argued that 
this compressibility decays with a universal power law as 
a Mott plateau is approached in a fermionic system. We 
investigate this quantity for a bosonic system. The value 
of I is somewhat arbitrary and was previously chosen to 
be larger than the correlation length in the insulating 
phase. Here we take the sum over the entire system, to 
remove the uncertainty associated with the size of I, and 
hereafter we remove the index I. The compressibility still 
behaves in essentially the same way as if we limit the sum 
to some smaller region. 

Extending the range of the summation makes the def- 
inition of the local compressibility identical to the def- 
inition used in Refs. jlOlll 1| . and the sum of the local 
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FIG. 6: Variance in the local-density Ai as a function of den- 
sity for the unconfined case [t/V = 0.05 (*), t/V = 0.1 (+), 
t/V — 0.2 (x)]. The solid line shows results for the confined 
system with a square potential. 



compressibility at all sites is the global compressibility. 
This compressibility diverges for the homogeneous sys- 
tem as the phase transition is approached from the su- 
perfluid phase, and is zero in the insulating phase, as 
mentioned above. Finite-size effects prevent the com- 
pressibility from diverging in a finite system, and instead 
the compressibility approaches zero as a smooth func- 
tion as the insulating state is entered. The way in which 
the local compressibility approaches zero in a finite but 
unconfined system is therefore a finite-size effect. 

In Fig. we show this decay of the local compressibil- 
ity as the Mott insulator is approached. It appears that 
the compressibility decays algebraically as we approach 
the Mott plateau, which is the behavior described for 
the fermionic system^. The exponent is very similar for 
the different strengths of the hopping. However, we also 
note that the compressibility for the homogeneous sys- 
tem follows the same behavior as the confined system, 
in contrast to the behavior depicted for the fermionic 
case. It appears that this algebraic decay is not the re- 
sult of universal behavior as a Mott plateau is approached 
in the confined system, but rather another property of 
the homogeneous model that carries over to the confined 
system. It is reasonable that the confined and the un- 
confined model display the same type of finite-size ef- 
fects, and therefore the observed similarity between the 
confined and the unconfined system is expected. In the 
bosonic case we find that the exponent that describes the 
decay has a numerical value of about 1.0, as compared to 
about 0.7-0.8 in the fermionic case£ The exponent ap- 
pears slightly smaller for the confined systems, but the 
statistical fluctuations are larger in this case, which prob- 
ably explains the small difference. The same exponent is 
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FIG. 7: Local compressibility K4 for the unconfined case 
[t/V = 0.05 (*), t/V = 0.1 (+), t/V = 0.2 (x)] and the 
confined system with a square potential [t/V = 0.05 (dashed 
line), t/V = 0.1 (solid line)]. 

obtained if we analyze the transition to zero density at 
the edge of the trap, which is expected since this is the 
same kind of phase transition. 

We have done a short calculation for the fermionic sys- 
tem and we found that for both the unconfined and the 
confined systems the local compressibility decays to zero 
in an algebraic manner. This is different from the pre- 
viously reported behavior™, described above, where only 
the confined system displays such a decay. A possible rea- 
son for the discrepancy is that we use a grand canonical 
ensemble, instead a of canonical ensemble, which seems 
to be more appropriate when observing particle number 
fluctuations. We found that the exponent for the un- 
confined model is very close to one and that the confined 
model may have the same, or a slightly, smaller exponent. 

By releasing the atoms in the trap the momentum dis- 
tribution, rifc, can be obtained experimentally^ This is 
probably the most accessible observable that can be stud- 
ied both experimentally and theoretically. The momen- 
tum distribution is the Fourier transform of the particle- 
particle correlation function, 

n k = J2exP iHri ~ ri H4 c i)- ( 5 ) 

In the homogeneous model the particle-particle correla- 
tion function is expected to decay algebraically in the 
superfluid phase, while it decays exponentially in the 
insulating phase. By a simple rescaling this has been 
demonstrated also for the trapped system^ 

In Fig. [5] we show the particle-particle correlation func- 
tion from a few selected points in the trapped system. 
We can clearly see the exponential decay of the correla- 
tion function within the insulating phase, and the much 
slower decay in the superfluid phase. We also note that 
correlation functions from points in the superfluid regions 
decay exponentially as they enter the Mott insulating re- 
gions. Especially interesting are points close to the phase 
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FIG. 8: The particle-particle correlation function (cjcj) for a 
few selected lattice points i to all other points j in the lattice. 
The dashed line represents the density profile, shown here 
to indicate the location of the Mott insulating phases. The 
dotted line marks the correlation function from a point on the 
phase separation boundary between the super fluid and Mott 
insulating regions. 



separation boundary. The dotted line represents the cor- 
relation function from such a point, which displays ex- 
ponential decay on one side and a long-range algebraic 
decay on the other side. 

To obtain the momentum distribution we first average 
the particle-particle correlation function over the whole 
system. In Fig. [5] the sum of all the particle-particle cor- 
relation functions is presented. This function is compared 
with the function obtained if the correlation functions for 
different ^-values for the homogeneous system are added. 
The chemical potentials are chosen with equal spacing in 
the whole range covered by the effective potential in the 
well. The sum can immediately be compared with the 
total correlation function for the system with a linear 
confining potential. However, for the other potentials 
the correlation functions have to be reweighted since we 
use equally spaced values of the chemical potential in the 
homogeneous system. To reweight the correlation func- 
tion we multiply the distribution with the derivative of 
the confining potential with respect to i, the lattice posi- 
tion. We notice that the agreement is very good for small 
distances, while the correlation function for the homoge- 
neous system has a longer tail, which can be expected, 
since the insulating phases cut off the correlations in the 
trapped system. 

When the correlation function is Fourier transformed 
to get the momentum distribution, the difference at large 
distances in the correlation function gives a difference 
in the momentum distributions for small k values. The 
momentum distribution for the confined systems is com- 
pared with the momentum distribution for the homoge- 
neous system in Fig. ^3 and except for the peak value 
we notice a reasonable overall agreement. 

Previous calculations have indicated that the appear- 




7 



3 





Quartic 




Square 




Linear 




Unconf. 







-150 -50 50 150 



12 r 

10 

8 




1 ■ — 1 ■ 1 ■ 1 ■ — 1 ■ 1 ■ — 1 

-3-2-10 1 2 3 

ka 



FIG. 9: The spatial average of the particle-particle correlation 
function for the confined system is compared with the value 
for the homogeneous system for t/V = 0.1. 
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FIG. 10: The momentum distribution for the confined system 
is compared with the value for the homogeneous system for 
t/V = 0.1. 

ance of a fine structure in the momentum distribution 
can be used to detect the formation of a Mott insulating 
plateau in the middle of the wel\m& In another study the 
fine structure is claimed to be a finite-size effect^ We ob- 
serve fine structure in the Fourier transform of correlation 
functions originating from individual sites in the super- 
fluid phase. This structure appears as a consequence of 
the abrupt cut in the correlation function obtained when 
the Mott insulating phases or the edges of the system 
are reached, as mentioned above. An example of this 
structure is shown in Fig. 1111 However, the sum of all 
the correlations is smooth since the cutoff in the correla- 
tion functions appears at different distances for different 
lattice points. In case a small lattice is studied some 
fine structure may remain and this supports the earlier 
conclusion^ that the satellite peaks are a finite-size effect. 



FIG. 11: An example of the fine structure when the particle- 
particle correlation function from one point in the superfluid 
phase is Fourier transformed. 



IV. CONCLUSIONS 

We examine the accuracy of a local-density approxima- 
tion for the one-dimensional Bose-Hubbard model using 
a quantum Monte Carlo method. In this approximation 
the behavior of the confined system is approximated with 
that of the unconfined system. In most of the interest- 
ing parameter regime one can obtain quantitatively accu- 
rate results using this approximation. In this manner the 
state diagram for a confined system can be determined 
to a high degree of accuracy from knowledge of the phase 
diagram of the homogeneous Bose-Hubbard model. We 
demonstrate how the approximation fails where the cor- 
relation length diverges in the homogeneous system, and 
also when the local density varies rapidly from site to site. 
Even for the case of the nonlocal momentum distribution 
the local-density approximation works fairly well. Con- 
trary to earlier studies we find no evidence of universal 
scaling as an insulating plateau is approached. Since the 
confined systems do not display a continuous phase tran- 
sition this is to be expected. Finally, we provide more 
evidence for the case that an observed fine structure in 
the momentum distribution is related to a finite-size ef- 
fect. 
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